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Three-body correlations for the ground-state decay of the lightest two-proton emitter ®Be are 
studied both theoretically and experimentally. Theoretical studies are performed in a three-body 
hyperspherical-harmonics cluster model. In the experimental studies, the ground state of ^Be was 
formed following the a decay of a ^"C beam inelastically excited through interactions with Be 
and C targets. Excellent agreement between theory and experiment is obtained demonstrating the 
existence of complicated correlation patterns which can elucidate the structure of "Be and, possibly, 
of the A=6 isobar. 

PACS numbers: 23.50-l-z, 23. 20. En, 21.60.Gx 



I. INTRODUCTION 

Two-proton {2p) radioactivity was predicted by V.I. 
Goldansky in 1960 fi| as an exclusively quantum- 
mechanical phenomenon. True three-body decay, in his 
terms, is a situation where the sequential emission of the 
particles is energetically prohibited and all the final-state 
fragments are emitted simultaneously. These energy con- 
ditions are illustrated for ^Be in Fig. [T] which shows that 
the '^Li ground state (g.s.) is not fully accessible for se- 
quential decay. Since the experimental discovery of "^Fe 
two-proton radioactivity in 2002 [1,13], this field has made 
fast progress. New cases of 2p radioactivity were found 
for ^"Zn ^, i^Mg [§], and, maybe, "*Ni M. The 2n cor- 
relations were recently measured in "^Fe 3, "'^^Mg [9l.[l0|- 
i^Ne [13, and i°C Very interesting [ll , but, so far, 
controversial [ll| case is possibility of 2p emission off de- 
formed isomeric state in ^"^Ag. All these decays exhibit 
complex correlation patterns. It is argued that studies 
of these patterns could provide important information 
about the structure of the decaying nuclei. 

With this active research as the background, there is 
one case which has been unduly forgotten. The ^Be nu- 
cleus is the lightest true two-proton emitter in the sense 
of Goldansky. As this is expected to be the simplest case 
(smallest Coulomb interaction, expressed cluster struc- 
ture with closed-shell core), a full understanding of its 
physics would provide a reliable basis for all future stud- 
ies of 2p decay. However until now, theoretical work on 
^Be was limited to predicting the energies and widths of 
its states. In addition, precise experimental data do not 
exist. The last experimental work dedicated to correla- 
tions in ^Be g.s. is about 20 years old. 

There is one more aspect which makes the ^Be case 
of special interest. In the last decade, large efforts have 
been directed to studies of ^He with special emphasis to 
the understanding of the halo properties in this compar- 
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FIG. 1: Energy levels and decay scheme for ^Be The ®Be 
g.s. is a true two-proton emitter in the sense of Goldansky: 
the sequential decay of this state is not possible as the lowest 
possible intermediate, ^Li g.s., is not energetically accessible. 



atively simple and accessible case. The associated litera- 
ture comprises hundreds of titles. To deduce the correla- 
tions in the neutron halo, one has to excite (e.g., Coulex) 
or destroy (e.g., knockout reactions) this nucleus. There- 
fore, the interpretation of the experimental data is influ- 
enced by the need to clarify details of the reaction mech- 
anism [47]. However, information about mirror system 
can be obtained without all this "violence" . The isobaric 
analogue state in ^Be decays to the a+p+p channel all 
by itself, providing the differential data on correlations. 
This data can be used directly to elucidate the structure 
of ^Be without the need to deal with the details of the 
reaction mechanism. Thus an important opportunity ex- 
ists for a better understanding of ^He properties through 
detailed studies of the ^Be. This has not been exploited 
previously. 

In this work, we provide detailed theoretical calcula- 
tions of the three-body decay characteristics of ®Be in a 
three-body cluster a+p+p model. We demonstrate that, 
in certain aspects, ^Be may be a preferable tool for stud- 
ies of the A=6 isobar, especially considering the high 
sensitivity of observables to the details of the theoreti- 
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cal models. We then discuss previous experimental and 
theoretical works on ^Be. Subsequently, we report on 
an experiment where ^Be fragments are formed after the 
Qf-decay of "'^^C projectiles excited by inelastic scattering. 
These data cover the complete kinematic space accessible 
for three-body decay and the correlations are compared 
to the theoretical predictions. 

The h = c = 1 system of units is used in this work. The 
following notations are used: Et is the system energy 
and Esr is the three-body resonance energy relative to 
the three-body a+p+p threshold. 



These Jacobi variables are given in "T" Jacobi system 
(see Fig. [2]). The hyperradial components x^^(p) of the 
WF equation[21 possessing the pure outgoing asymptotics 

^^+\p,n,) = £ x^+l {>cp) j^^ (n,), (5) 

are matched to approximate boundary conditions of the 
three-body Coulomb problem obtained in Ref. [l^. The 
radial components of this WF at large p values can be 
represented as 



II. THEORETICAL MODEL 

The theoretical framework of this paper is largely the 
same as that developed for the two-proton radioactivity 
and three-body decay studies in Refs. [13, [H, [H, 113, US, 
[2^ . It is based on the three-body cluster model using 
the hyperspherical-harmonics (HH) method. The pre- 
dictions obtained with this approach were found to be 
in very good agreement with experimental widths and 
momentum distributions 

In this section, we sketch the necessary formalism 
emphasizing only the points which differ from previous 
treatments. 



A. Hyperspherical harmonics method 

For narrow states, the time-dependent wavefunction 
(WF) in a finite domain can be parameterized as 



(1) 



The radial part of this WF can be obtained with good 
precision as a solution of the inhomogeneous system of 
equations 



(6) 



In general, the functions Tif)^}^ are some linear combi- 
nations of Coulomb functions with the outgoing asymp- 
totic G + iF. The functions Jk^^ {^p) are hyperspheri- 
cal harmonics coupled with spin functions to total spin 
J. "Multyindex" 7 denote the complete set of quan- 
tum numbers except the principal quantum number K: 
7 = {L,S,lx,ly\- The value /Cmax truncates the hyper- 
spherical expansion. The hypermoment x is expressed 
via the energies of the subsystems E^, Ey or via the Ja- 
cobi momenta k^, ky conjugate to Jacobi coordinates X, 
Y: 
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(7a) 
(7b) 



= 2MEt = 2M{Ex + Ey) ^2kl + -kl, (7c) 
= {Ok, ,^ky} , Ok^ aictanlEx /Ey] . (7d) 



A more detailed picture of the "T" and "Y" Jacobi sys- 
tems in coordinate and momentum spaces can be found 
in Fig. H 

The set of coupled equations for the functions x'"^'' has 
the form 



H-E3r) ^-L^^ (p, ^p) = -I (r/2) *box(p, ^p) , (2a) 



(2b) 



Ef — T -\- Vcp {^crti ) ^cp {^cp2 ) Vpp {^PlP2 ) 



Here ^'box and E^r are the eigenfunction and the eigen- 
value of the equation 



(3) 



^f-^ *box(p,f^p)-0. 



solved with a "box" boundary condition at large p. The 
hyperspherical coordinates are defined via the Jacobi vec- 
tors 
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{9p, fij., ily} , dp = arctan 




(4c) 
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2M ^ Vk^^k'y {p)x^kK ip) + * TM Xk^ (p) , (8) 

K'Y 

where C ~ K + 3/2 is "effective angular momentum" and 
VK'j.K'-y'ip) is "three-body potential" (matrix elements 
of the pairwise potentials); 



dQp (Qp) J2 V^, (r., ) (np) (9) 



and 



*box(p,f^p) = P~'^'J2^K''^p')^k''^^p^ ■ (It^) 
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FIG. 2: Independent "T" and "Y" Jacobi systems for the 
core+A^+A'^ three-body system in coordinate and momentum 
spaces. 



B. Width and momentum distribution 

Equation [2] is first solved with an arbitrary value of 
r and then the width is found according to the "natu- 
ral" definition as the flux j through a hypersphere with 
large radius Pmax divided by the internal normalization 
N ("number of particles" inside the sphere): 
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if7 



dn, 
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dp 



The differential flux through the hypersphere pn 
fined as 

djipma.x,flp) _ 



(11) 
(12) 

(13) 
is de- 
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(14) 



If, for sufficiently large p, the coefficients Af^^^f" in Eq.[6] 
become independent of p, then the coordinate distribu- 
tion becomes identical to the momentum distribution, 
i.e., 



dfln 



dfl. 



(15) 



Further discussions of the validity of this approximation 
(Eq. I15|) , and detailed expressions for the momentum dis- 
tributions, can be found in Ref. Il6l. 



C. Potentials 

The NN potential is taken either as a simple s-wave 
single-Gaussian form BJ (from the book of Brown and 
Jackson (Hi) 



with Vo = —31 MeV and rp = 1.8 fm, or the realis- 
tic "soft-core" potential GPT (Gogny-Pires-de Tourreil 

M)- 

The Coulomb potential of the homogeneously charged 
sphere rgph = 1.852 fm is used in the a-p channel. In 
addition for this channel, we use an ^-dependent potential 
SBB (Sack-Biedenharn-Breit ^S]) 

V{r) = I/W exp(-rVr2) + (£ • s) Ve, exp(-rVrg) , (17) 



system where tq = 2.30 fm, yi"^ = 50 MeV, Vc 



.(1) 



Vi^^ = -23 MeV, and V/, = -11.71 MeV. Historically, a 
somewhat modified SBBM potential has been used in the 
calculations of A=6 isobars in order to better reproduce 
the binding energies (e.g., Ref. [l^). Later it was reafized 
that it is more consistent to provide the phenomenolog- 
ical binding-energy correction using an additional short- 
range three-body potential (see, e.g., the discussion in 
Ref. [llj). In this work, we used a short-range three- 
body potential of the form 

Vsip) = Sk-,j<'yV^''^ /[I + exp((p - po)/d3)] , (18) 

where po — 2.5 fm and da = 0.4 fm. This "short-range" 
three-body potential (note the small difluseness) does not 
distort the interactions in the subbarrier region which 
was found to be important for consistent studies of decay 
properties. 

Three sets of nuclear potential are employed in 
this work. They are denoted as PI (SBB+BJ), P2 
(SBB+GPT), and P3 (SBBM+GPT). The values ofvj°^ 
used with potential sets PI, P2, and P3 are —11.14 MeV, 
-13.22 MeV, and -0.64 MeV, respectively. Throughout 
this paper when the potential set is not specified, the 
results of the calculations with the P2 set are shown. 



-47.32 MeV, 



V{r)=Voexp{-r'/r'„) 



(16) 



D. Reaction models 

In general, different definitions for the width of a de- 
caying state coincide only in the limit when the width is 
very small. For the ground state of ^Be, the definition of 
Eq. [TT]is not very precise, as this state is comparatively 
broad (P = 92 ± 6 keV) and thus the internal normaliza- 
tion N (Eq. I13p is sensitive to the integration limit pi^t . 
For reasonable values of pmt ranging from 10 — 20 fm, 
the uncertainty in the width (Pnat) is about 25% [see 
Sec. nil CI Fig. [TWa)]. This problem does not exist for 

narrow 2p emitters (P < 1 eV) where the WFs x'k-I ^'"^ 
vanishingly small under the Coulomb barrier. The den- 
sities for the dominating components of the ^He and ^Be 
WFs are shown in Fig. [31 For ^Be, it is clear that the 
WF under the barrier is not negligible. 

For moderately broad states, there are alternative ways 
to derive the width. These involve either the study of the 
3^3 scattering or the study of a particular reaction. For 
technical reasons, the latter is preferable for our applica- 
tion. For example, in order to determine the population 
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FIG. 3: Densities IXK-yip)]^ and |xi?7(p)P for the largest com- 
ponents of the ®He and ^Be g.s. WFs. 



In an adiabatic approximation, we can assume that the 
radial part of kinetic energy is small under the centrifu- 
gal barrier in the channels where this barrier is large and 
can be approximated as a constant. In this approxima- 
tion, the FR procedure is reduced to the construction of 
effective three-body interactions by the matrix 

operations 

{K + 3/2){K + 5/2) 



Ef - E 



5Kj,K'-i' (22) 



where 



of ^Be in a charge-exchange reaction on ^Li at zero angle, 
Eq. [2] can be reformulated as 

(h-Et) ^iil{p,%) = 

E-rE-f'^*^'Lf(^'^p)- (19) 

i M 

This notation is based on the fact that for angles close to 
zero, the transitions in charge- exchange reactions, in the 
limit of high energies, are provided by the Gamow- Teller 
operator. Although this reaction is different to the one 
studied experimentally in this work, it is sufficient for our 
computational purposes. Namely, we will demonstrate 
that for the ®Be g.s. population, the choice of the reaction 
mechanism is not very important (there are still some 
exclusive situations, which wc will discuss elsewhere). 

Using the source function of Eq. the cross section 
for the population of the three-body continuum is pro- 
portional to the outgoing flux of the three particles on a 
hyperspherc of some large radius p = Pmax'- 

da{ET)/dET ^ j{pm.^, ^p) ■ (20) 

Differentials of this flux on the hyperspherc provide an- 
gular and energy distributions among the decay products 
at the given decay energy Et in analogy with Egs. 1141 and 

E. "Feshbach" reduction 



Summations over indexes with the bar are made for the 
eliminated channels. We typically eliminate the channels 
with K > KpR, where Kpn provides the sector of the 
hyperspherical basis where the calculations remains fully 
dynamical. We take the "Feshbach energy" Ef in our 
calculations as Ef = E. 

There are two ways to control the reliability of the 
FR procedure, (i) The "soft" method is to vary KpR 
from the maximum attainable in the dynamic calcula- 
tions downwards for fixed /^max- The results, in princi- 
ple, should coincide, (ii) The "safe" method is to take 
Kmay. in the range attainable for dynamic calculations 
and compare the "reduced" i^max — * KpR calculations 
(with much smaller dynamic basis size Kpn) with com- 
pletely dynamic calculations with i^^max- For ^Be, these 
considerations show that we can safely use KpR = 14. 
However, the even safer value of KpR — 22 is used in 
this work. 



III. GROUND STATE 

There are several convergence characteristics that 
should be understood before reliable results on ^Be are 
obtained. The convergence character is quite different for 
all the observables of interest and also depend strongly 
on the interaction in the p-p channel. 



Although the HH calculations for ^Be can be per- 
formed with -f^max = 22 — 26, these basis sizes may not be 
sufficient to obtain good convergence for all observables. 
However, the basis size can be effectively increased using 
the adiabatic procedure based on the so-called Feshbach 
reduction (FR) [2l| . Feshbach reduction eliminates from 
the total WF "if = '^p + '^g^ an arbitrary subspace q using 
the Green's function of this subspace: 

E[p = Tp + Vp ~ VpqGqVpq . (21) 



A. Convergence of energy and width 

Because of the problem mentioned in Sect. Ill Dl we 
need to begin our studies with the energy dependence 
of the cross section. The convergence of the cross-section 
profile with increasing size of the basis is demonstrated in 
Fig. m The main character of the convergence is clearly 
seen here; the centroid energy decreases, while the width 
grows significantly. 
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FIG. 4: (Color online) Energy profile of the ®Be g.s. pop- 
ulated in the charge-exchange reaction with ®Li. The re- 
sults are shown as a function of the basis size Jfmax where 
Kfr ~ 22. For ifmax < AVji and no Feshbach reduction is 
needed. 



FIG. 5: (Color online) Dependence of "Be g.s. width on the 
decay energy Et- Predications are shown for the three po- 
tential set P1-P3. The dotted curves show the approximation 
of Eq. [251 



The cross section for the ^Be g.s. population, shown in 
Fig. m clearly has a profile close to a slightly asymmet- 
ric Lorentzian. Can the profile of this three-body reso- 
nance be described by appropriately-modified R-matrix 
type expressions? A curious result is obtained here, the 
resonance profile, shown in Fig. [4] by the solid curve, can 
be fit with amazing precision by the following expression: 



TiEj 



<j{Et) 



To 



T{Et) 



Et\^ + (1 - a) f ^ 



(24) 



(25) 



where Fo = 98 kcV and a = 0.65. Equation [Ml is the or- 
dinary expression for the inelastic cross section of an iso- 
lated resonance. The parameterization of Eg. 1251 was cho- 
sen because, for the single-channel penetration through 
the hyperspherical barrier with K = 0, the energy depen- 
dence of the width can be inferred as T(Et) ^ E^. For 
if = 2 one has T{Et) ~ E^ (see, e.g., Ref. It should 
be understood that the K = component is equivalent 
to a "phase volume" with the characteristic energy be- 
havior of ^ E"^). The energy dependence of the width 
obtained by Eq. [21] almost coincides with the calculated 
dependence of this width in a reasonable energy range 
(see Fig. [5] when one uses a=0.63 and 0.52 for potential 
P2 and P3, respectively). If we take the actual calcu- 
lated partial widths for the K — and 2 components 
from Table |Vl then the value of a can be estimated as 



Nk=q/{Nk=o + Nk=2) « 0.58 



(26) 



This is quite close to the value 0.65 obtained by a fit. 

The existence of this simple approximation, despite the 
fact that there are Coulomb interactions and other nu- 
merous channels involved, may demonstrate that the dy- 
namics of the ^Be g.s. decay is largely defined by the 
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FIG. 6: (Color online) Convergency of the resonance energy 
Eir and the width F for the "Be g.s. as a function of the basis 
size Kuia.^ • 



penetration through the hyperspherical barriers. Pos- 
sibly, this is due to the comparatively large ^Be decay 
energy of E'^r — 1.371 MeV. Simple estimates shows that 
the state is "sitting" somewhere straight on the top of 
the Coulomb barrier. 

It was found that the value of j(i?T) for ^Be g.s. is not 
sensitive to the particular choice of the source in Eq. [191 
which is typically within the width of the line [48] . This 
means that the width defined by the procedure of Eqs.[Ml 
and [2S1 is very reliable. We can fine tune the value pmt 
in Eq. [131 so that the definition of the width in Eq. [Til 
coincides with the definition in Eq. [Ml and subsequently 
we can reliably use Eq.[Tll All of the potential sets P1-P3 
needed pint ~ 12.5 fm. 

The convergences of the predicted resonance energy 
and width as a function of the hyperspherical basis size 
are shown in Fig. [HI for each of the potential sets. In 
all cases, our calculations are fully converged. The res- 
onance energies are forced to approach the experimental 
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FIG. 7: (a) Sensitivity of the width as defined in Eg. [TTlto 
the size of the "internal region" pint ■ The dot shows the value 
of Pint at which this width coincides with that defined via 
the cross-section profile Eqs. [191 [IHl [24] and [25] (b) Relative 
precision of the width as a function of the matching radius 

Pmax • 

value E^r = 1.371 MeV. This is achieved by fine tuning 
the phenomenological potential of Eq. [181 this is a neces- 
sary approach in order to provide reasonable predictions 
for the decay characteristics. We can see that while the 
calculations with PI and P2 (SBB potential in the a- 
p channel) arc in good agreement with each other and 
with the experimental value, the width obtained with P3 
(SBBM potential) is far too large. 

An expected feature observed here is the much slower 
convergence of the calculations with a realistic potential 
in the NN channel. An important, but often disregarded 
fact, which one can see in Fig. [S] is the much slower con- 
vergence of the width as compared to the energy. This 
means that, in general, an energy convergence does not 
guarantee the convergence of other important character- 
istics. As we will see in Sect. IIII C[ the situation with 
momentum distributions is even more complicated than 
it is for the widths. 

The sensitivity of the width to a number of the other 
parameters in the calculations is demonstrated in Fig. [T] 
Figure [7fa) shows the sensitivity of the width defined by 
Eq. [TT] to the size pint of the region where the internal 
normalization is calculated. The stability of the calcu- 
lations to the dynamical range Pmax is demonstrated in 
Fig. [TJJb). To attain 1% numerical precision in the width 
calculations, we need to go beyond 60 fm in the hyperra- 
dius p. 



B. Features of the momentum distributions in Be 

The correlations in the decay of ®Be include both the 
generic features of the 2p decays, as discussed earlier in 
Refs. Bain, and some peculiarities which we present 
in more detail now. For nuclear states with J < 1/2 (as is 
the case for ^Be g.s. decay), the three-body correlations 
can be completely described by 2 parameters. There are a 
total of 9 degrees of freedom for three particles in the final 
state. Of these, three describe the center-of-mass motion, 
three describe the Euler rotation of the decay plane (for 




FIG. 8: (Color online) Complete correlation picture for ®Be 
g.s. decay, presented in "T" and "Y" Jacobi systems. 

J < 1 all its orientations are quantum-mechanically iden- 
tical), and the three-body decay energy is fixed. Thus we 
are left with two parameters to describe the correlations. 
It is convenient to choose the energy distribution param- 
eter £ between any two of the particles and the angle 9k 
between the Jacobi momenta: 

e = E,/ET , cosiek)^^f^. (27) 

These parameters can be constructed in any Jacobi sys- 
tem and for ^Be there are two "irreducible" Jacobi sys- 
tems, called "T" and "Y" , see Fig. [H The distributions 
constructed in different Jacobi systems are just different 
representations of the same physical picture. However, 
different aspects of the correlations may be better re- 
vealed in a particular Jacobi system. 

Predictions for the complete correlation picture of ^Be 
g.s. decay are shown in Fig. [Sj for both the "T" and "Y" 
Jacobi systems. Schematic figures are included in this fig- 
ure to help in visualizing the correlations associated with 
different regions of the Jacobi plots. The main features 
of these distributions are: 

1. The energy distribution in the "T" system has a 
double-humped profile which is an indication of 
the [p^] configuration dominance which was pointed 
out in very early papers on ^Be [l^, [2^, Ho]- This 
double-humped configuration is expressed more in 
coordinate space (see the internal region in Fig.[T51) 
and only marginally "survives" in the asymptotic 
region. The internal peaks in Fig. [T3] have the 
special names of "diproton" (protons are close to 
each other) and "cigar" (protons are in-line with 
a-particle) configurations [26|. 

2. There are kinematical regions where the presence of 
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FIG. 9: (Color online) Dependence of energy distribution be- 
tween the proton ("T" system) in the decay of ^Be g.s. on 
the decay energy Et- 



particles is suppressed due to Coulomb repulsions. 
Strong suppression in the a-p channel in regions 
(b) and (d) and a smaller suppression in the p-p 
channel in region (e) are predicted. 

3. There are enhancements due to the p-p final-state 
interaction in regions (a) and (f). The ^Li g.s. res- 
onance in the a-p channel is not accessible for de- 
cay. However, some hint of its presence can be ob- 
tained from the enhancement in region (g) . This is 
a "back-to-back" configuration, where protons fly 
in the opposite directions. However, the reason for 
the enhancement of such a conflguration is not fully 
understood. 

4. The angular dependence in the "T" system al- 
most vanishes for regions (a) and (c) {E^/Et ^ 
and Ex/Et '-^ 1). It is clear that in the limit 
Ex/Et — > and E^/Et — > 1 the dependence on 
the relative orientation of kj; and kj, should be- 
come degenerate. However at intermediate values 
of Ex/Et, this dependence is very pronounced. 



The total-energy distribution in the "Y" system 
(see Fig. [TU] for the projected distributions) is al- 
most a symmetric bell-shape. This is the energy 
distribution between the core and one of the pro- 
tons and its symmetry reflects the symmetry be- 
tween protons. In heavy two-proton emitters, this 
distribution becomes very narrow and almost com- 
pletely symmetric. 
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FIG. 10: Sensitivity of the energy and angular distributions 
in the decay of ®Be g.s. to the choice of the potential set. 
Results are shown for both the "T" and "Y" Jacobi systems. 



should be suppressed at low energies. The pure [s^] con- 
figuration should produce a featureless "phase-volume" 
energy distribution 



dj/dEx - y/ExiEr-Ex) 



(28) 



The sensitivity of the projected distributions to the 
choice of the potential set P1-P3 is demonstrated in 
Fig. [ini The angular distribution in the "T" system and 
the energy distribution in "Y" systems are practically in- 
sensitive to this choice. The other projected distributions 
demonstrate sensitivity on the level of 10 — 15%. How- 
ever, local differences in certain kinematical regions are 
much larger. 

Figures [51 and [TUl demonstrate what we call the "soft- 
ness" of the ^Be system: minor variations in the condi- 
tions or computational details lead to a noticeable vari- 
ations in the observable properties. Heavier 2p emitters 
appear to be much "stiffer" in this respect. 



The correlation predictions shown in Fig. [5] are ob- 
tained on resonance. The dependence the of energy cor- 
relation on the decay energy of ^Be is demonstrated in 
Fig. O The double-humped shape of this spectrum be- 
comes less pronounced when the energy decreases. With 
smaller energy, the relative contribution of the [s^] con- 
figuration to the decay grows compared to the [p^] con- 
figuration. The latter has an additional centrifugal com- 
ponent to the barrier and its contribution to the width 



C. Convergence of the momentum distributions 

In our calculations there are two projected distribu- 
tions which are practically insensitive to convergence is- 
sues (the angular distribution in the "T" system and the 
energy distribution in the "Y" system). The other two 
distributions (the angular distribution in the "Y" system 
and the energy distribution in the "T" system) demon- 
strate strong sensitivity. The convergence of the energy 
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FIG. 11: (Color online) Convergency of the "T" energy dis- 
tribution in the decay of ®Be g.s. as a function of the basis 
size Tfm ax ■ 



distributions are illustrated in Figs. fTT] and fT2l 

The convergence of the energy distribution between 
protons has a very curious character, see Fig. [11] From 
^max = 8 to i^max — 22 this distribution is very sta- 
ble [several curves almost coincide, see Fig. [Qja)]. Then 
from -ftTjjiax — 24 to A'niax ^ 70 the distribution changes 
qualitatively, and up to -fCmax ~ 100 there is still a notice- 
able variation [Fig. [Hl^b)]. Hopefully with JC,„ax = 110, 
we have a well converged distribution. Calculations with 
small basis sizes (e.g., -ftTmax < 70) for ^Be should pro- 
vide a qualitatively wrong energy distribution in the "T" 
system. Similarly for the angular distribution in the "Y" 
system. 

This "softness" of the ^Be system makes it a very com- 
plicated object to study. Minimum basis sizes which 
provide convergence for the energy and width are far 
from sufficient for calculations of momentum distribu- 
tions. This is a feature which we probably do not face 
in heavier 2p emitters as the Coulomb interaction in the 
core-p channel plays a more dominant role in the decay 
dynamics. 

The radial convergence of the energy distributions is 
illustrated in Fig. [TH Calculations with pmax < 300 fm 
are clearly insufficient to stabilize the distribution. How- 
ever by Pmax = 800 fm, the distributions seem to be 
well converged. Could there be some noticeable modi- 
fications of the distributions due to further propagation 
in the long-range Coulomb field? This question was an- 
alyzed in Ref. [l6| for ''^Fe using the classical trajectory 
approach. The complete stabilization takes place in "^^Fe 
at /9 ~ (3 — 6) X 10** fm, with a major part of the effect 
originating at p < 1 x 10^ fm. The decay energies of ^Be 
and "^^Fe g.s. are similar and the Coulomb interaction is 




EJEr 

FIG. 12: (Color online) Dependence of energy distribution in 
the decay of ^Be g.s. on the maximal dynamic range of the 
calculation pmax- For the "T" Jacobi system, Ex is energy 
between two protons and in the "Y" Jacobi system, is 
energy between core and one of the protons. 

~ 12 times weaker in ®Be. Therefore, the majority of the 
long-range effects should be taken into account in calcu- 
lations with pmax ~ 1000 fm. The ^Be calculations of 
this work were typically done with pmax = 1200 fm. 



D. Structure of the '^He and '^Be g.s. 

From another point of view, one can benefit from the 
"softness" of ^Be system. The high sensitivity of the ob- 
servables to the details of the model ingredients increase 
our ability to discriminate these features and hence im- 
prove our ability to elucidate the details of the nuclear 
structure. 

Detailed information about the ^He and ®Be g.s. WFs 
is provided in Table IVl In general, there is high degree of 
isobaric symmetry between the ^He and ^Be WFs in the 
internal region. This is not true, however, for the K = Q 
component, which differs the most. The reason for this 
is shown in Fig. [3] where the magnitude of the K = Q WF 
in asymptotic region is comparable to its magnitude in 
the nuclear interior. Hence the nuclear boundary is not 
defined for this component in ^Be. This is also seen in 
Table [H which provides the information about the WF 



TABLE I: Weights of the shell-model-like configurations {f'] 
in the ^He and ®Be g.s. WFs in percent for the Jacobi "Y" 
system. The normalizations of the ®Be components are found 
for integration radius pint = 12.5 fm. 













f^Bc 






PI 


P2 


P3 


PI 


P2 


P3 




8.11 


8.58 


8.35 


10.54 


11.15 


10.84 




90.91 


90.30 


90.37 


87.98 


87.18 


87.17 




0.47 


0.53 


0.61 


0.69 


0.77 


0.95 




0.41 


0.43 


0.50 


0.60 


0.65 


0.77 
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FIG. 13: Spatial correlation densities F)p for the ®He and ^Be g.s. WFs in the "T" system. For ''Be, two panels provide 

the view in different radial ranges. Pay atention to difference in the scales. The variation from the top to the thick contour 
line in the ®He panel corresponds to the whole scale variation in the two ®Be panels. 



in approximate "shell model" terms. After looking at 
the radial behavior of the WF's components in Fig. [3l 
we find that the concept of isobaric symmetry is relevant 
here strictly speaking only for the most interior region 
of the WF (p < 4 — 5 fm) . Beyond this point the radial 
behavior in ^He and ^Be differ drastically. 

The weights of the components in Tables |V] and |T] are 
in very good relative agreement for the different poten- 
tial sets P1-P3. Evidently these major features of the 
structure are not that sensitive to the fine details of the 
interactions. 

It can be seen that the partial widths of the ^Be 
WF components in Table |V] are drastically different as 
compared to weights Ni in the internal region. This is a 
reflection of complicated dynamics in 2p decays, the WFs 
are strongly "rearranged" in the subbarrier region and by 



TABLE II: Radial properties of the He g.s. WF and some 
observables obtained for ®He and ®Be g.s. with potentials PI, 
P2, P3. 



value 


PI 


P2 


P3 


Exp. 


(P> (fm) 


5.088 


5.156 


5.491 




(rjvjv) (fm) 


4.482 


4.502 


4.884 




(rcjv) (fm) 


4.113 


4.172 


4.430 




{tm) (fm) 


3.211 


3.248 


3.469 




(r.) (fm) 


1.321 


1.171 


1.232 




r-mat (fm) 


2.396 


2.421 


2.540 


2.30 ± 0.07 [35] 










2.48 ± 0.03 [36] 


rch" (fm) 


2.103 


2.012 


2.048 


2.054 ±0.014 [33] 


rch' (fm) 


2.113 


2.043 


2.079 


2.068 ±0.011 [34] 


SGT^He^^Li) 


5.004 


5.058 


4.930 


4.745 ± 0.009 [2] 


AScoui (MeV) 


2.351 


2.302 


2.111 


2.344 [2] 


rC^Beg.,.) (keV) 


98 


112 


154 


92 ± 6 [2] 



"Theoretical values in this row are obtainGd using the generally 
accepted value for the neutron charge radius ^'^^{n) = —0.1161 fm^. 

''Theoretical values in this row are obtained using r^Y^{n) = 0.012 
fm2 [si. 



the long-range Coulomb pairwise fields. The ^He and ^Be 
WF correlation densities are shown in Fig. [13] The WFs 
are nearly identical in the internal region, while in the 
asymptotic region for ^Be we can clearly see how this 
"rearrangement" is taking place. Comparing different 
potential sets P1-P3 in Table El we see that PI and P2 
calculations are almost identical, while the major partial 
width in P3 differs strongly. We conclude that the decay 
dynamics is mainly defined by core-p interaction. 

Geometric properties of the ^He g.s. WF and several 
observables obtained for ^He and ^Be g.s. are shown in 
Table [iTl The root mean square values are given for p, 
tnn (distance between valence nucleons), Tcat (distance 
between nucleon and core) , (distance between valence 
nucleon and cm.), (distance between core and cm.). 
The differences between these geometric characteristics 
for PI and P2 are typically around 1%. In the case of 
P3, the differences are significantly larger. The Bgt val- 
ues obtained with P1-P3 also agree within 1.5%, but all 
differ more from the experimental value. Here, the "ex- 
perimental" Bgt value is obtained using the ^He lifetime 
Ti/2 = 806.7 ± 1.5 ms 0, and the /3-decay constants of 
/<(0+ ^ 0+) = 3072.40 s and A = 1.268. It has al- 
ready been discussed in the literature that the 4 — 7% 
disagreement here could be connected with both the WF 
quality and the renormalization of the weak constant 26] . 
Therefore, we give no definite conclusion about quality 
of the models here. 

The next most precisely known characteristic for ^He 
is its charge radius. Recent studies have defined Vch with 
increasing precision (ssl . [3^ ] . The relative uncertainty of 
this value is now about 0.5%, while variations in the cal- 
culated value are around 4% for P1-P3. However, com- 
parison of this value with those theoretically calculated 
is not completely model independent. The theoretically 
calculated charge radius of ^He is noticeably sensitive to 
the neutron charge radius. The latter is inferred theoret- 
ically, rather than measured experimentally. This means 
that there exists considerable systematic uncertainty in 
the determination of the charge radii. According to our 
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estimates, this uncertainty can be as large as 2%. This 
fact somewhat relieves the constrains on the WF con- 
nected with this observable. One can see in Table Hi] 
that the P2 and P3 calculations, containing realistic NN 
potentials can be regarded as consistent with the exper- 
iment. 

The matter radius of ^He is defined in the cluster 
model using the matter radius of the a-particle. The 
value rniat(a) = 1.464 fm is derived from the charge ra- 
dius Teh (a) = 1.671 using the neutron and proton charge 
radu; r^in) = -0.1161 fm^, rch(p) = 0.875 fm. The ex- 
perimental data on matter radii have large a systematic 
uncertainty. This is probably the reason for the contro- 
versial signal obtained in different experiments (see two 
examples in Table . This observable so far does not 
seem to have discriminative power for theoretical models. 

The Coulomb shift Ai^coui and ^Be g.s. width obtained 
with PI and P2 are in a good agreement with experiment. 
Some overestimation of the width in the three-body clus- 
ter model can be expected due to the admixture of dif- 
ferent configurations in ^Be WF. The weight of such ad- 
mixtures can be estimated as 6 — 14%, based on the PI 
and P2 widths. However, the Coulomb shift and width 
obtained with P3 are clearly not acceptable. Our over- 
all feeling is that the cumulative information on ^He and 
^Be g.s. is sufficient to choose P2 as the only acceptable 
potential. 



IV. THEORETICAL DISCUSSION 

As we have already mentioned, most of the attention in 
the studies of the A=6 isobar has been paid to ^He. Even 
in the studies of ^Be, there are only few works which stud- 
ied it's width. In addition, there has been are only limited 
studies of the ®Be g.s. decay correlations. The first con- 
sistent calculations of the ^Be three-body decay width 
were performed in Ref. [32| using the integral formalism. 
In papers p^ . [l5| . the quantum-mechanical formalism 
for two-proton radioactivity and Coulombic three-body 
decay studies was developed. In these papers, the inte- 
gral formalism was criticized in application to the decays 
of systems with strong three-body Coulomb interactions 
and a more preferable way to calculate widths was pro- 
posed [see, Eq. [11] . The value F = 90 keV was obtained 
in Ref. [ij] with the PI potential (i^max — 20), which as 
we can see in Fig. [SI is reasonably well converged. 

In our approach, the effects of antisymmetrization are 
taken into account in a simplified way. However, there are 
studies that treated the ^Be decay as a 6-body problem. 
In RGM calculations [111, the ^Be width of F = 160 keV 
for E^r = 1.52 MeV was found using the complex scaling 
method. Scaling this value to the experimental 2p decay 
energy with the help of Fig. [5| we obtain F = 125 keV 
which is considerably larger than the experimental value. 
An interesting algebraic method was developed for stud- 
ies of ^Be decay in Ref. Here, the hyperspherical 
decomposition is used for the WF both in the internal 



region (6-body HHs) and in the asymptotic region (three- 
body cluster HHs) . A calculated width of F = 72 keV was 
obtained for = 1.172 MeV which scales to F = 110 
keV at the experimental 2p decay energy. In addition, 
we can expect a 10 — 15% reduction due to the absence 
of the S — I component in these calculations. This com- 
ponent is important in the internal region, but does not 
contribute to the width significantly. In addition, we 
can also expect roughly a factor of 2 increase due to the 
small basis size (i^max = 10) used in the asymptotic re- 
gion in Ref. [13] ■ According to Fig. [HI with A'^ax = 10 
we can ex pec t only 60% of the width, at most. It seem 
that Ref. [40| is more a concept demonstration, rather 
than a realistic calculation. Therefore at the present mo- 
ment, it is not possible to draw any conclusions about 
importance of the 6-body effects in calculations of the 
^Be decay properties. 

The width of the ^Be g.s. was calculated in Ref. fssj via 
a method analogous to ours (hyperspherical harmonics) , 
but having certain technical diff'erences. An approximate 
treatment of the 3 — > 3 scattering is introduced in this 
work and the width is extracted from the energy behavior 
of the phase shifts. The width obtained was F = 65 keV 
for Esr = 1.26 MeV which scales to F = 84 keV at the ex- 
perimental 2p decay energy. It can be found in Ref. jsl] 
that the calculation does not seem converged. If we ex- 
trapolate from F = 84 keV using the convergence curves 
for P2, P3, then the value F = 110 keV is obtained, which 
is is a good agreement with our P2 result. 

An important result of the present work is the clear 
demonstration that any approach purporting to give sat- 
isfactory description of the ^Be g.s. decay properties 
should have a certain "dynamic range" both in radial and 
functional spaces (see Table IIIIl) . It can be found that not 
all of these conditions are satisfied in these other works 
dedicated to ^Be. 

Our calculations demonstrate a noticeable sensitivity 
of the observables in the decay of ^Be g.s. to the ingredi- 
ents of the model. Table [TTj demonstrates that this sen- 
sitivity is enhanced in ^Be compared to ^He. Typical 
variations of the observables for ®Hc are 0.5 — 4%, while 
in ^Be there is about a 60% difference in between the 
widths calculated with PI and P3. The tunneling process 
can be seen as a kind of a "quantum amplifier" , which 
drastically emphasizes minor features in the structure. 
For that reason, it is possible that the indirect probe of 
^Be decay is a more sensitive tool for determining the 
halo properties of ^He than direct investigations of ^He 
itself. We are referring to precision measurements of the 
correlations in ^Be decay which are discriminative with 
respect to the fine details of the momentum distributions. 
In the experimental studies presented in this work, the 
quality of the data is approaching fulfilment of such a 
high precision request. 
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FIG. 14: Experimental energy distributions between protons 
in the decay of ^Be measured in (a) Ref. [12] and (b) Ref. [30| . 
The theoretical prediction (PI) is provided only to guide the 
eye, as now we have no idea about the required experimental 
corrections. 



V. EXISTING EXPERIMENTAL KNOWLEDGE 
ABOUT '^'BE 

Very precise results about the energy and width of the 
^Be g.s. were obtained in the early studies: Et — 1371(5) 
keV, r = 89(6) keV [41]. The current value of the width 
is only slightly different T = 92(6) keV 0. 

The first measurements of ^Be decay correlations were 
made in Ref. [12], see Fig. [UT a). They determined the 
energy spectrum of a-particles reconstructed in the ^Be 
cm. frame. For ^Be g.s. events, this spectrum is the 
same as the correlation spectrum between two protons. 
The authors could not fit the data using simplistic decay 
scenarios (phase volume, diproton decay, simultaneous 
emission of p-wave protons) and concluded: "...no inco- 
herent sum of the processes considered here will fit the 
data. Perhaps a full three-body computation is necessary 
to understand the energy spectrum." 

This ground-state decay, as well as decays of the 2+, 
T = states of the A^Q isobar, was further investigated 



TABLE III: Minimal dynamical ranges of calculations re- 
quired to provide reasonably converged different observables 
for ^Be. Different basis sizes are required for simplistic BJ 
and realistic GPT potentials in the p-p channel. 



value 




r 


distributions 


Pmax (fm) 


20 


60 


300 


i^max (SBB+BJ) 


16 


30 


80 


/STmax (SBB-FGPT) 


40 


70 


110 



in the series of works by the Kurchatov Institute group 
[H, [m, and Refs. therein], see Fig.[T4tb). They de- 
veloped a method of analyzing the p-p correlations in 
the framework of a three-body partial-wave decomposi- 
tion and applied this to the three-body decays of light 
nuclei (29l.[30l]. In particular, the first kinematically com- 
plete study of ^Be proved the existence of three-particle 
p+p +a correlations with S{p-p) = 1 and S(p-p) = 
[SOllsH which matched the three-body components found 
theoretically in the p-shell structure of ^Be ^2^. One of 
the important result for ^Be g.s. was the realization that 
S{p-p) — and S{p-p) — 1 components of the WF should 
produce very different correlation patterns. The presence 
of an "admixture" of S{p-p) ~ 1 component to the WF 
was demonstrated by an experiment performed with spe- 
cial kinematics. In these works, the concept of "demo- 
cratic decay" was coined. This describes the specific de- 
cay mode for three-body systems, when the events are 
not highly focused in narrow kinematical regions, but are 
distributed broadly ("democracy" among different kine- 
matical regions). "Democratic decay" is now a popular 
term for this class of phenomena, but the correlations in 
^Be decay have never been studied since that time. The 
spectra shown in Figs. [14] (a) and (b) are not in complete 
agreement with each other. Furthermore, there are large 
statistical uncertainties and the geometry of experiments 
may cause cuts in kinematical space which make compar- 
ison the theory difficult. It is clear a modern experiment 
on ^Be decay was needed. 



VI. EXPERIMENT 
A. Experimental Method 

The Texas A&M University K500 cyclotron facility was 
used to produce a 200 pnA beam of ^"B at E/A — 15.0 
MeV. This primary beam impinged on a hydrogen gas cell 
held at a pressure of 2 atmospheres and kept at liquid- 
nitrogen temperature. A secondary beam of E/A = 
10.7 MeV was produced through the i''B(p, n)i'^C 
reaction and separated from other reaction products us- 
ing the MARS spectrometer This secondary beam, 
with intensity of 2 x 10^ s~^, purity of 99.5%, an energy 
spread of 3%, and a spot size of 3.5 x 3.5 mm was in- 
elastically excited due to interactions with 14.1 mg/cm^ 
Be and 13.4 mg/cm^ C targets. Ground-state ^Be frag- 
ments were created from the a decay of these excited 
particles. Following the decay of the ^Be g.s. fragment, 
the final exit channel is 2p-|-2a. 

The four decay products were detected in an array of 
four Si E-AE telescopes located in a plane 14 cm down- 
stream of the target. The telescopes, part of the HIRA 
array Q , consisted of a 65 /.tm thick, single-sided Si-strip 
AE detector followed by a 1.5 mm thick, double-sided Si 
strip E detector. All Si detectors were 6.4 x 6.4 cm in area 
with their position-sensitive faces divided into 32 strips. 
The telescopes were positioned in a square arrangement 
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FIG. 15: Contour plot showing the distribution of the two 
possible ^Be excitation energies that can be associated with 
the two 2p+Qf subevents ordered by the maximum and min- 
imum values. The dashed line indicates the threshold for 
which correct identification of the 2p-\-a subevent associated 
with "Be decay is achieved in the simulations. The ridge as- 
sociated with ®Be g.s. decay is indicated by the arrow. 



FIG. 16: (Color online) Experimental distribution of Et for 
selected events is shown by the data points. The solid curve 
indicates the distribution predicted by the Monte Carlo sim- 
ulation with the addition of a smooth background (dashed 
curve). The G6 gate used to select "Be g.s. events and the 
two gates (Gsi, Gb2), used to estimate the background under 
the peak, are indicated. 



with each telescope offset from its neighbor to produce 
a small, central, square hole through which the unscat- 
tered beam passed. With this arrangement, the angular 
range from = 1.3 to 7.7° was covered. More details of 
the experimental arrange can be found in Ref. [ll[. 



B. Monte Carlo Simulations 

Monte Carlo simulations of the experiment were per- 
formed in order to determine the experiment bias and 
to understand the effects of the gates applied to remove 
unwanted 2p+2a events. The simulations included the a 
decay of the parent fragments and the correlations 
between the ^Be decay products are sampled according 
to the theory of Sect. IIII Bl The effects of energy loss 
and small-angle scattering of all the decay products were 
considered following Refs. [isLlioj. 

Simulated events were passed through a detector fil- 
ter and the effects of the position and energy resolution 
of the detector were added. The "detected" simulated 
events were subsequently analyzed in the same manner as 
the experimental data. The velocity, excitation-energy, 
and angular distributions of the parent ^''C states were 
chosen such that the secondary distributions that passed 
the detector filter were consistent with the experimental 
results. Similar simulations for other decay modes were 
found to reproduce the experimental resolution [ll| . 



C. Event Selection 

Apart from a-^Be g.s. decay, these are many other ^"C 
decay modes that lead to the 2p+2a exit channel and 
thus the detected events must be suitably gated to re- 




E*("'C) [MeV] 



FIG. 17: (Color online) Experimental distribution of ^''C ex- 
citation energy for a-^Be events selected in this study. 



move these unwanted decays. Of particular importance 
is the rejection of the large yield of decays where the ^"C 
fragments undergoes two-proton decay (either sequential 
through ^B or prompt) leading to the creation of an ®Be 
g.s. [ll|. These events can readily be identified from the 
correlations between the two a particles. The distribu- 
tion of relative energy {E"^) between the two a particles 
contains a strong, narrow peak corresponding to ^Be g.s. 
decay fu\. This peak has a FWHM of 38 keV and sits 
on a negligible background [ll| thus allowing for a clean 
rejection of these events with the gate E^"^ < 0.2 MeV. 
Our Monte Carlo simulations suggests this gate has es- 
sentially no significant effect on true a-^Be g.s. decays 
with only 0.01% of detected events being rejected. 

The remaining events have contributions from a-^Be 
g.s. andp-^B {E*=2A3 MeV) decays The latter ^B 
excited state does not decay through ^Be g.s. but under- 
goes a three-body decay like the ^Be ground state. For 
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both of these decays modes, there is a difficultly is try- 
ing to find the intermediate state (either ^Be or ^B) as 
there are two possible ways to construct this fragment 
from the detected 2p+2a exit channel. Let us concen- 
trate on the ^Be g.s. fragments first where we must de- 
termine which of the two detected a particles was the 
one initially emitted from the ^"^C parent and which was 
produced in the decay of ^Be. To this end, the ^Be 
excitation energy for the two ways of constructing the 
^Be fragment are determined and ordered according to 
their maximum and minimum values; i?*(^Be)iiiax and 
i?*(^Be)ijiin. A two dimensional plot of these two exci- 
tation energies is shown in Fig. 1151 A prominent ridge 
centered around i?*(^Be)min — corresponding to ^Be 
g.s. decay is clearly visible. For those events in this ridge 
structure, the identification of which a particles was pro- 
duced in ®Be decay is clearly the one associated with 
£^*(^Be)„,in when E*{^Be)n,^y, > E*{^Be)n,in- However 
when £'*(^Be)max ^ -E*(^Be)„iin the Monte Carlo simu- 
lations indicate that misidentifications will occur. These 
simulations suggests that for i?*(^Be)niax — ^^*(^Be)min = 
0.5 MeV, the probability of misidentifying the a particles 
is 0.03%. This condition is indicated in Fig. [15] by the 
dashed line and only events above this line were used in 
the subsequent analysis of the experimental data. One 
can see from Fig. [15] that this condition does not signifi- 
cantly cut into the ridge structure and the Monte Carlo 
simulations suggests we lose 4.7% of the remaining a-^Be 
g.s. events with this gate. 

The remaining ridge structure still sits on a back- 
ground. Part of this background can be traced to ^°C^ 
p+^B{E* — 2.43 MeV) decays. These events can be iden- 
tified from -E*(^B)max and E*{^B)„^in information in a 
manner similar to the a-^Be g.s. events. A ridge struc- 
ture also is evident in this case and it also sits on an 
non-negligible background, which in turn has contribu- 
tions from a-^Be g.s. decay. Although one cannot com- 
pletely separate all p-^B and a-^Be events, we do reject 
events in the i<^*(^B)niin ridge structure. This results in 
a slightly diminished yield of true a-^Be g.s. events, but 
more importantly, it reduces the relative background un- 
der the ^Be ridge structure shown in Fig.[T2] The Monte 
Carlo simulations suggests only 2.7% of the remaining 
true a-^Be g.s. events were rejected by this condition. 

The distribution of Et for the final selection of events 
is shown in Fig. \W\ by the data points. The FWHM 
width of the peak associated with ^Bc g.s. is 220 keV 
which is larger than the intrinsic value of F = 92 keV 
due to detector resolution. The solid curve indicates the 
simulated distribution after a smooth background con- 
tribution (dashed curve) is added. This simulated dis- 
tribution reproduces the experimental results quite well 
confirming that the Monte Carlo simulations correctly 
model the experimental resolution. Figure [TBI also shows 
the gate G6 used to select ^Be g.s. fragments and the 
two gates, Gbi and Gb2 which, when combined, were 
used to estimate the background in the G6 gate. In all 
subsequent results, this background has been subtracted. 




FIG. 18: Comparison of (a,c) experimental (exp) and (b,d) 
theoretical correlations for ®Be g.s. decay presented in both 
the (a,b) "T" and (c,d) "Y" Jacobi systems. The theoreti- 
cal distributions include the effects of the detector bias and 
resolution. 



The excitation-energy distribution of ^"C fragments as- 
sociated with the selected events is shown in Fig. [T7] 
There is localized strength around E*(^°C) = 7 MeV and 
a continuous distribution up to approximately 15 MeV. 
Thus many ^°C excited states are contributing to the 
detected ^Be g.s. yield. 



VII. COMPARISON OF THEORY AND 
EXPERIMENT 

Comparisons of experimental and predicted correla- 
tions in both the "T" and "Y" Jacobi systems are shown 
in Fig. [TH] The experimental results [Figs. [TSTa) and 
[TSf c'l] has been background subtracted and, for the pre- 
dicted distributions [Figs. [T5f b) andlTSTd)]. the effects of 
the detector resolution and bias has been incorporated 
via the Monte Carlo simulations. In this and subsequent 
plots, the simulated results has been normalized to the 
same number of counts as for the experiment data. In 
determining the Jacobi coordinates, there are two ways 
of choosing the order of the proton. For the experimen- 
tal events, Jacobi coordinates were determined for both 
of these ways and thus each event contributes two counts 
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FIG. 19: (Color online) Comparison of the experimental (data 
points) and predicted (curves) distributions of Ex/Et in the 
"T" (c), (d) and "Y" (a), (b) Jacobi systems. The blue dashed 
curves show the primary predicted distributions while the red 
solid curves include the effect of the detector bias and res- 
olution. The theoretical results were obtained with the P2 
potential. 
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FIG. 20: (Color online) Comparison of the experimental 
(data points) and different predicted (curves) distributions 
of S^/Bt in the "T" (c), (d) and "Y" (a), (b) Jacobi sys- 
tems. The red, blue, and green curves correspond to PI, P2, 
and P3 potential sets respectively. The effect of the detector 
bias and resolution is included for the theoretical curves. 



to the spectra. For "T" system, this forces the cos(6'fc) 
distribution to be symmetrized around cos {9k)=0. Gen- 
eral overall agreement between theory and experiment is 
found, although statistical fluctuations are the limiting 
factor for the experimental data. 

To allow for a more detailed comparison, we compare 
projections of the correlations on both the E^/Et and 
cos {9k) axes in Fig [121 The experimental data are in- 
dicated by the data points while the dashed and solid 
curves show the predictions before and after the sim- 
ulated bias of the experimental apparatus is included. 
Interestingly, the "soft" observables (energy distribution 
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FIG. 21: (Color online) Comparison of experimental (data 
points) and predicted (curves) cos{9k) distributions in the "T" 
(left) and "Y" (right) Jacobi systems for the indicated gates 
on E^/Et parameter. The red, blue, and green curves cor- 
respond to PI, P2, and P3 potential sets respectively. The 
effect of the detector bias and resolution is included. 



in the "T" system and the angular distribution in the 
"Y" system) which have the most sensitivity to the in- 
gredients of the theoretical calculations and its numerical 
implementation also have the largest bias induced by the 
detector apparatus. The other projected distributions 
(angular distribution in "T" and energy distribution in 
"Y" ) are practically unaffected the detector response. 

The same comparison of theory and data for all three 
potentials P1-P3 is shown in Fig. [201 All three sets of 
predictions reproduce the experimental data reasonable 
well. To highlight more details of the correlations, we 
show the cos [Ok) distributions gated on three equal re- 
gion of E,^/Et in Fig. [H] for the "T" and "Y" Jacobi 
systems. Reasonable agreement between the experiment 
(data points) and the three calculations (curves) is also 
found, although the PI and P2 calculation are somewhat 
better. To quantify this, we determine the per de- 
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TABLE IV: P^r degree of freedom for fits to the complete 
correlations data in the "T" and "Y" system with the three 
assumed potentials 



potential 




"Y" 


PI 


1.29 


1.25 


P2 


1.17 


1.14 


P3 


1.58 


1.45 



gree of freedom {x^ /v) of the theoretical fit to the two- 
dimension data of Fig. [THl These values are listed in 
Table HV] for both the "T" and "Y" systems. For a good 
fit were need /v ^1 and clearly both PI and P2 satisfy 
this criteria. Again we find the P3 calculation is some- 
what worse. 



VIII. CONCLUSIONS 

The first detailed studies of the correlations from the 
decay of ^Be g.s. are performed both experimentally and 
theoretically. We have found that certain correlations 
(namely, energy correlation between two protons and an- 
gular correlations in "Y" Jacobi system) are quite sen- 
sitive to the details of structure and interactions. We 
demonstrated that relative sensitivity of correlation pat- 
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TABLE V: Weights Ni of the dominating components of the He and Be g.s. WFs and the partial widths of the Be g.s. WF 
in percent. The results are for the Jacobi "T" system. The normalizations of the ®Be components are found for integration 

radius pint = 12.5 fin. 
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